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Braneworld modifications to the Friedmann expansion law can have an important effect 
on the cosmological evolution of the early universe. In particular, the primordial particle 
■ abundances crucially depend on the rate at which the universe expanded at early times. In 

(N ■ 

^ ' this article, we study the production of stable and unstable gravitinos, both from thermal 

I creation and from the decay of a heavy scalar, in the Randall-Sundrum II braneworld context. 

, We conclude that, depending on the value of the 5D fundamental Planck mass, some of the 

04 ■ usual standard cosmology constraints on the reheating temperature and on the mass of the 

heavy scalar can be evaded. 

§ ; PACS numbers: 98.80.-k, 98.80.Cq, 04.50. +h, 04.65. +e 
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^ . I. INTRODUCTION 

Q^' The (over)production of gravitinos can have striking impact on cosmology [1]. If gravitinos are 

(D ! 

r~| , unstable, their overabundance would spoil the success of big-bang nucleosynthesis (BBN). On the 



X 
J3 



other hand, if they are stable particles, which is the case if they are the lightest supersymmetric 
particle (LSP) contributing to the total amount of dark matter, the present abundance of gravitinos 
is constrained in order to avoid the overclosure of the universe. 

It is long known that gravitinos can be thermally produced through scatterings in the plasma. 
Howeve. .t has been .ecently pointed out Q.B. in .he context of the n„dnB p.ob.en, ||. that the 
decay rate of the scalar field into a pair of gravitinos is much higher than previously thought [5|] . 
This "new" source of gravitinos coming from a scalar field decay might have important consequences 
for the evolution of the universe. In fact, we expect that the universe was once dominated by a 
scalar field, for example, the inflaton (see, e.g., Ref. \^ for a review), or the dilaton and moduli 
fields in superstring theories. The production of gravitinos from infiaton/dilaton/scalar decays, 
and the problems that might result, have been widely studied in standard cosmology (SC) 
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In the last years, there has also been much interest in the so-called braneworld cosmology 
(BC) scenarios. Whilst theories formulated in extra dimensions have been around since the early 
twentieth century, recent developments in string theory have opened up the possibility that our 
universe could be a 1+3-surface, the brane, embedded in a higher-dimensional space-time, called 
the bulk, with all (minimal supersymmetric) standard model particles and fields trapped on the 
brane, while gravity is free to access the bulk (see Ref. [13t | for a review). A remarkable feature 
of BC is the modification of the expansion rate of the universe before the BBN era^. In the so- 
called Randall-Sundrum II (RSII) braneworld construction [3], in which one has a single brane 
with positive tension A, the Friedmann equation receives an additional term quadratic in the 



density 
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22 



23l |. Setting the 4D cosmological constant to zero and assuming 



that inflation rapidly makes any dark radiation term negligible, one obtains 



Sir 
3M 



Mp 



3 M| 



(1) 



where Mp = 1.22 x 10^^ GeV is the 4D Planck mass and M5 the 5D fundamental mass. Notice 
that Eq. ([T]) reduces to the usual Friedmann equation, H oc y^, at sufficiently low energies, /o <C A, 
but at very high energies one has H cc p. Successful BBN requires that the change in the expansion 
rate due to the new terms in the Friedmann equation be sufficiently small at scales ~ O(MeV); 
this implies M5 > 40 TeV. A more stringent bound, > 10^ TeV, is obtained by requiring the 
theory to reduce to Newtonian gravity on scales larger than 1 mm. 



As is well known, the above modification in t 
implications for inflation [24], baryogenesis 
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le Friedmann expansion rate can have important 



duction 



26l |. dark matter [2^ and thermal gravitino pro- 



23]. In this paper we aim to investigate its effect in the decay of a heavy scalar field into 



a pair of gravitinos. After a brief discussion on the decay of heavy scalars and the production of 
gravitinos in the RSII brane scenario, we study the resulting constraints on the parameter space, 
both for the stable and unstable gravitino cases. 



II. HEAVY SCALAR DECAY AND RSII COSMOLOGY 



We consider here the case in which the scalar particle X dominates the energy of the universe 
when it decays, p ~ px- This scenario can be achieved if the scalar is, for example, the inflaton. 



^ There are also brane constructions in which the expansion rate is modified after the BBN and gravity is modified 
at large distances, as in the Dvali-Gabadadze-Porrati model [14], which is a possible explanation for the present 
accelerated expansion of the universe. 
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We will assume that X is a singlet under the standard model gauge group. 

Of particular interest is the reheating temperature, T^h, which is defined here by assuming 
an instantaneous conversion of X energy into radiation, when the decay width of X, Tx, equals 
the expansion rate of the universe. If the total energy density of the scalar X is instantaneously 
converted into radiation, then we can identify px = PR = (vr^/30) , where 5* is the effective 
number of relativistic degrees of freedom at the temperature T. In particular, g^, = 915/4 = 
228.75 in the minimal supersymmetric standard model (MSSM), for temperatures above the SUSY 
breaking scale ~ O(TeV). For a radiation dominated universe, Eq. ([T|) can then be written as 

4" 



where 



-'I 



(2) 



is the usual SC expansion law, and we have defined 



8-lxlO^GeV , (4) 



— ^ = ».i X iU' LieV — -7— — 

as the transition temperature from the high energy brane regime expansion to the standard ex- 
pansion^. From the condition 

H{Trh) = Tx (5) 

one gets the reheating temperature Trh, where Tx is the total decay rate of X. In the high-energy 
limit of BC, i.e. for T > , Eq. Q reads 

27r35, 



45 M| 



(6) 



and the condition in Eq. (0) leads then to the relation 

1/4 



In the SC limit (T <C Tt) one obtains 

/ 45 \l/4 



^ To compute the numerical value, we have used the value of as in the MSSM. Hereafter, the same value will be 
used to obtain other numerical estimates in the present work. 
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The total decay rate of X is determined from all the interactions of X with the other particles 
of the supersymmetric model, and so it is model dependent. However we can in general write 

rx = a,^. (9) 

In fact, when X has only Planck suppressed interaction, at is expected to be of order one. For 
example, if X is the heavy moduli field that couples to gauge supermultiplets in the gauge kinetic 
function through a Planck suppressed interaction, the decay rate of X is of the form of Eq. ([9]) 
with at = 0(1) 

Using Eq. Q, the reheating temperature in the BC regime can be written as 

Trh = 6.8 X 10 GeVa/ ( ^T^TnTTTTT ) ( TTTTnT^ ) ' (^0) 
while in the SC limit one obtains 



IQiOGeVy VlQiOGeV 



Mx \ 



3/2 



r,, = 5.7xlO^GeVaW-,^ . (11) 



IQio GeVy 

At this point one should notice that T^-h > 7 MeV in order to be consistent with BBN theory. 
This implies that in BC regime the relation M^Ms > 4.8 x 10^° '^^^ GeV^ should be satisfied, 

2/3 

while in SC one should have Mx ^ 2.5 x 10 a^ GeV. However, in order to have the BC regime 
one should also impose T^h > Tt, i.e. to obtain such a low reheating temperature in the BC 
regime one needs also a low Tt which corresponds to a low value of Ms. Such low values of Ms 
are not compatible with the bound Ms > 10^ TeV, obtained by requiring the theory to reduce 
to Newtonian gravity on scales larger than 1 mm. This implies that the only relevant bound is 
coming from the SC regime 

Mx > 2.5 X 10^ GeV . (12) 
One should also point out that in order to have Trh > Tt, one must impose 

Mv>(^) M.. (13) 

Here we are mainly interested in the gravitino production, and hence in the decay of a heavy 
scalar field X into a pair of gravitinos 

X 1p3/2 + i^3/2 ■ (14) 

Before we proceed with the discussion of the decay of the heavy scalar into gravitinos, we 
should specify better our brane setup, namely the localization of gravitinos on the brane. We 



are considering here the supersymmetric extension of the RSII braneworld model [29] . As already 
mentioned, in the Randall-Sundrum construction under consideration, only gravity propagates in 
the bulk: the zero- mode graviton is localized around the ultraviolet (UV) brane (our brane), while 
the graviton Kaluza-Klein (KK) modes are localized around the infrared brane (which in the RSII 
is located at the infinite boundary). In the supersymmetric generalization, all the fields in the 
same supermultiplet will take the same configuration in the fifth dimensional direction. Hence, 
the zero-mode gravitino will also be localized around the UV brane. As an approximation we will 
consider it as a field residing on the UV brane, as done previously in other works 



2$. The 



zero-mode gravitino will obtain its mass through the super-Higgs mechanism, as in the usual 4D 
supergravity. Since supersymmetry is broken only locally on the UV brane, the KK gravitinos will 
remain the same. We will not consider here the production of these KK gravitinos and its possible 
influence on the brane particle content. We will only consider the zero-mode gravitino as a field 
residing on the brane and apply 4D supersymmetry and the usual Boltzmann equation to describe 
its production. 

Let us now briefly discuss the decay channel in Eq. (1140 . For simplicity, we will assume that 
Mx » ^T^3/2• In fact, the gravitino is likely to be much lighter than either a heavy modulus^ or 
the inflaton, because a too large gravitino mass requires a flne-tuning in the Higgs sector due to 
the anomaly-mediated effects. With this assumption, the decay rate for this process, as recently 
~\, can be written as 



calculated in Refs. 



"3/2 

= ^Ml- ^^^^ 
The coupling constant 03/2 is defined by the relation 

with G being the total Kahler potential, G = K + hi\W\^ , where K and W are the Kahler potential 
and the superpotential, respectively. The subscript X {X) denotes differentiation with respect to 
the X (X) field and ( ... ) stands for the VEV. In general one has Q, E, Q, 13, l3, U] 

where the value of the constant c depends on the explicit model. For example, in the scenario 
discussed in Ref. one has c = \/3- For a moduli field, 03/2 is expected to be order unity Q], 



^ Many proposed solutions to solve the moduli problem require the moduli to have large masses, e.g. Mx > 
100 TeV t30|. 



6 



and hence the decay of the moduh field into a pair of gravitinos can be significant, thus leading to 
cosmological disasters. If X is the inflaton, then 03/2 can take values much smaller than unity, as 
in some new and hybrid inflationary models. If 03/2 is zero then no gravitinos are produced in the 
decay of X. This is what happens in the simplest chaotic inflationary models, for which Gx = 
in the vacuum (see Ref. 0, 9, 10|). In the following we will leave as a free parameter. 

One should note that the expression in Eq. (|15p for the decay rate of the gravitino cannot be 
applicable for H > 771,3/2. As pointed out in [8], the decay proceeds only if the Hubble parameter 
is smaller than the gravitino mass, since the chirality flip of the gravitino forbids the decay to 
proceed otherwise. This imposes a bound on the maximum value that Mx can take, and also on 
the maximum reheating temperature. The condition H{Trh) < Tn3/2 is, by the definition of T^/j 
[cf. Eq. ([5])], the same as Tx < Hence 

2.5 X 10^ GeV af^' < Mx < 5.3 x lO^^ af^' {jp^) ' (18) 



where the first inequality is the one already presented in Eq. (jl2p . The upper bound on the 
reheating temperature turns out to be 

in the high-energy brane regime, and 



r.,<2.2xl0i°GeV(^^^)'^' (20) 



in the SC limit. 

Summarizing, one would expect that the brane effects will change the gravitino production for 
reheating temperatures in the range 



or, equivalently, heavy scalar masses in the interval 

2_\^/^ _ „ -5/3 / "^3/2 \l/3 



{£) M, < Mx < 5.3 X 10^3 af^' {j^) • (22) 



III. GRAVITINO ABUNDANCE 



As referred in the Introduction, a viable cosmological supergravity-inspired scenario has to avoid 
the so-called gravitino problem. For unstable gravitinos, their decay products should not alter the 
BBN predictions for the abundance of light elements in the universe. On the other hand, if they 
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are stable particles they will contribute to the dark matter content, and thus one has to avoid the 
overclosure of the universe. 

In the following we will estimate the contributions to the gravitino abundance, both from 
thermal production and from the decay of a heavy scalar. Then we will discuss the constraints 
coming for stable or unstable gravitinos. 



A. Brane cosmology vs standard cosmology 

Gravitinos can be produced by the decay of the heavy scalar particle X, and, during the 
reheating phase, they are also thermally produced through scatterings in the plasma. The total 
abundance is then given by both contributions 

^3/2 = = Y3/2 + ^3/2 ' (23) 

where s{T) = 2iT'^g^T^ / A5 is the entropy density, and stands for the contribution from 

thermal scatterings {X decay) to the total gravitino abundance^. 

Let us estimate the contribution of thermal scatterings at reheating. The gravitino abundance 
at a given temperature T < Trh is obtained from the Boltzmann equation 



'■3/2 



dt 



+ 3i/n3/2 = C3/2(T) , (24) 



where C^/2{'^) is the collision term given by [34] 



C3/2(T) ^ a{T) 1 + PiT)^ —J , (25) 



"^3/2/ 

where nig is the low-energy gluino mass; a{T) and /3(T) are slowly-varying functions of the tem- 
perature. Assuming constant entropy, the integration of Eq. (I24p yields 

^2 



(26) 



with 

3/2 ^,(n^\ ( rtii 



/ 45 \ a(T) / mf, \ 

^(T;™3,„™,)=(-) Jd.{l + ,3(T)^J ■ (27) 



In fact, there are other potential sources of gravitinos in addition to these two The gravitinos can be produced 
in the decay X — > X + ^p3^2, where X is the fermionic partner of X, however for this to happen is required a 
large mass hierarchy between X and X [3l|. Also, the decay of X usually produces superparticles, followed by 
cascade decays into gravitinos. Finally, gravitinos can be produced before the decay of X in such a way that their 
abundance remains sizable after the dilution by the reheating. In this work, we will not consider any of these 
contributions. 



8 



This is an approximate solution, obtained by integrating the Boltzmann equation in the BC regime 
for Trh > T > Tt and the SC Umit for T < Tj, considering a small dependence of A{T;m^i2-,fng) 



on the temperature. In Ref. 



23], an analytic solution in terms of a Gauss hypergeometric function 



was found. However the result presented in Eqs. (j26j) and (p7|) is, for our purposes, a very good 
approximation . 

In the high energy limit of BC, Eq. (I26p leads to 

Yl/l^^2A{Tum^l2,m~g)Tt , (28) 

while in the SC limit one obtains the usual result 

^3/2 - ATrh-, , m-g) Trh ■ (29) 

In the BC regime, the abundance of gravitinos produced by thermal scatterings is proportional 
to the transition temperature instead of the reheating temperature. This has been proposed as 



a possible so^ 
leptogenesis 



ution for the thermal gravitino problem in the context of supersymmetric thermal 
in SC, the gravitino bound imposed on the reheating temperature, Trh S 10^ — 



10^'^ GeV (depending on the gravitino mass and the hadronic branching ratio 8|, lid . l33l|). is 
barely compatible with the simplest models of supersymmetric thermal leptogenesis which need 
Trh > 10^ - 10^ GeV. 

Let us now look at the contribution of the decay process ()14p to the total gravitino abundance 
[cf. Eq. ([23]) ] . The amount of gravitinos produced by this process can be written as 

n% = 5%.§^. (30) 



where 



r_o /9 



is the branching ratio of the decay channel (|14p . If 03/2 ^ 1 (or /and at ^ 1), -B3/2 can be much 
smaller than unity. 

The abundance of gravitinos produced from the decay of X can then be written as 

_ «3/2 Trh _ 4/2M'xTrh 



Notice that, in the brane scenario, H{Trh) is in fact a function of Trh and T^ (or M5). In the BC 
and SC limits, this equation leads, respectively, to 



Mx Yf Mg Yf Trh \ 



y3%- 8.9x10-1^42 ^^010 Gevj ^0^° GeVy* ^lO^GiVy' 

^-3/2«y^(^)'V^)"'\ (33) 
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a3,,M,(GeV) 



FIG. 1: (Color online) Contours of constant ¥3^2 = 10^^^, 10^^^, lO^^'"^, 10^^"' (from lighter to darker green) 
in the Trh — ct3/2^Ix and q;3/2 — Mx planes, for BC with — 10^^ GeV (full lines) and SC (dashed lines). 
The (black) horizontal dashed line in the left panel represents the value of Tt; the (black) vertical line in the 

1 /3 

right panel shows the value of Mx — {2 /of) M5. The (blue) dash-dotted line represents the lower bound 
on Trh from -B3/2 < 1 for ^3/2 = 10^^". The gravitino mass is taken to be 7713/2 — 1 TeV and, in the right 
panel, we fixed at = 1. The maximum values for T^h and Mx shown in the plots are the ones given by the 
upper bounds of Eqs. (pij) and (P^ . respectively. 



and 



X 



3/2 



1.4 X 10"** 



Mx V ( Trh y' 



8.6 X 10-6 S3/2 at 



Mx 



1010 GeV 



1/2 



(34) 



Note that in both regimes the gravitino abundance is inversely proportional to Trh, in opposition 
to the thermal abundance in SC. 

In Fig. [T]we present the contours of constant in the Trh — "3/2-^x and 03/2 — Mx planes. 
We show the SC case (dashed lines) and the BC one (full lines), for a fixed value of M5 (M5 = 
10^1 GeV). We have chosen the gravitino mass to be 7713/2 = 1 TeV. In the right panel, we see the 
change of behavior in the contours for Mx > (2/aj) M5, for which the high-energy regime of 
BC is valid. This regime occurs for Trh larger than Tt, as we can see on the left panel. For the 
contour plot in the 03/2 — Mx plane we have chosen at = 1. One should notice that if we fix at at 
a smaller value, then the lines of constant would move downwards, i.e. only smaller values of 
Q^3/2 would be allowed. 

In the left panel plot, the lower bound on Trh from -B3/2 < 1 is also shown by choosing 03/2 = 
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FIG. 2: (Color online) Same as in Fig. [H but for the total abundance of gravitinos: ¥3/2 = 
10^^^, 10^^^, 10^^^, 10^^** (from lighter to darker green). In the computation of the thermal abundance, 
the gluino mass was taken to be m-g = 1 TeV. 



10 (dash-dotted line). This bound can be derived from Eq. (|3T]) . which can also be rewritten as 



Br, 



a 



3/2 M'x 



'mH{Trh)Ml ' 
and that in the high-energy hmit of BC reads 

Ms \V Mx 



^3/2 = 5.9xl0-^a^/,, GeV 



1010 GeVy VlO^GeVy 



This leads to 



T,/, >4.9a 



1/4 ( Ms V'^ f a3/2Mx 
3/2 VlQiOGeVy V GeV 



3/4 



In SC one has 



D nil n-5 2 f \ ( Tj-h 

53/, = 9.1x10 o.,,J-^\ {^^^ 



and hence 



Trh > 3.0 X 10"** a 



1/2 / gg/aMx y/^ 
3/2 I 105 GeV j 



(35) 



(36) 



(37) 



(38) 



(39) 



Of course, if one requires a smaller value for -63/21 ^ larger value of the reheating temperature is 
needed. And if 03/2 is larger, the bound will be less stringent, in both cases, for a fixed 03/2 Mx- 
In Fig. [2] we show the contours of constant 1^3/2 = ~^ ^3/2 same plane as in Fig. [U 

assuming rrig = 1 TeV in the computation of the thermal contribution to the gravitino abundance 
(hereafter we will consider this mass for the gluino). In the right panel, we see that in the SC 
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1 /3 

regime, i.e. for Mx < {2/al) ' M5, the heavy scalar mass Mx has an upper bound (coming from 
the thermal production contribution), which corresponds to the upper bound on the reheating 
temperature (see left panel). These bounds can be easily derived from Eq. ([29]) . which can be 



rewritten as 



^ 10- 



16 



10-22 GeV 



^1 



106 GeV 



5.7 X 10"^*^ at 



10-22 GeV-1 



Mx \ 
1010 GeV J 



3/2 



(40) 



from which one gets 



Trh < 10^ GeV 



A{Trh;rn3/2,mg) 
10-22 GeV-^ 



1 /yBBN' 



3/2 



10-16 



(41) 



and 



Mx < 6.7 X lO^'' GeV 



N-l -2/3 /vBBN\ 2/3 



10-22 GeV 



-1 



10-16 



(42) 



On the other hand, in order to avoid the overproduction of gravitinos by X decays, the coupling 
constant 03/2 has to satisfy 



a3/2<2.1 X 10-5 a^' 



Mx 



lOio GeV 



or, in terms of the branching ratio. 



B^/2 < 1.2 X 10 



ari 



Mx 



-1/4 /vBBN\ 1/2 
f ^3/2 \ 

10-16 



-1/2 /vBBN' 
/ -^3/2 



(43) 



(44) 



* VlQiOGeVy \^ 10-16 
as already derived in Qj. 

1 /3 

For Mx > (2/0?)'/^ M5, or equivalently Trh > Tt, the high -energy regime of BC holds. In 
this regime, the upper bound on the Trh disappears. Indeed, this bound is due to the thermal 
production of gravitinos, and, as can be seen from Eq. (j28|) . the abundance of gravitinos produced 
thermally is proportional to Tt, 



yl% ^ 2 X 10 



-16 



~ 1.6 X 10 



A{Tt]m^l2,mg) 
10-22 GeV-1 
_17 \ A{Tt;m^i2,mg 
10-22 GeV-1 



106 GeV 

M5 



lOio GeV 



3/2 



(45) 



For m-g = 1 TeV, mg/a = 1 - 10^ TeV and T = 10'* - lO^^ GeV one has ^(T; mg/a, mg) ~ lO'^^ GeV'^ 
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Hence, in the high-energy regime the bound on the reheating temperature is replaced by a bound 
on the transition temperature, 

'A{Tt; 1713/2, rrigy^ 



Tt<5x 10^ GeV 

10-22 QeV" 

or, equivalently, on the 5D fundamental mass M5 

'A{Tt;my2,mg) 



' vBBN ■ 
-^3/2 

10-16 



(46) 



M5 < 3.4 X 10^° GeV 



n -2/3 /vBBN\ 2/3 
/ -^3/2 \ 



10 



-16 



(47) 



10-22 GeV-i 

This is also reflected in the disappearance of the upper bound on the mass of the heavy scalar 
imposed by the limits on the gravitino abundance, and the only bound that remains on this quantity 
is the one in Eq. (jlSp . which comes from requiring m3/2 > H. In fact, in this regime, in order to 
avoid overproduction of gravitinos in the X decay, one must have 



J,, ^4.5x10 ^^o5GeVy' ^10^° GeV j 



' vBBN ~ 
^3/ 2 

10 



-1/3 



-16 



(48) 



easily obtained from Eq. (|33|) . This means that, like in SC case, the X decay into gravitinos 
imposes a lower bound on the reheating temperature. In SC one must have [cf. Eq. ^4 



/ \ 2 /vBBN\ -1 

> 1.4 X 10^ f ^^^"l 
- '-^^ ' 105 GeVy ' ' 



10 



-16 



(49) 



The lower bound coming from -B3/2) presented in Eqs. p7p and (j39p . can be stronger than those in 
Eqs. (|18]) and (|19]) . depending on the value of B3/2 and/or 03/2- One finds that the bounds from 
-B3/2 can be disregarded, if one considers sufficiently large Q!3/2- 

We should also notice that in the high energy regime the upper bound on the coupling constant 
Q^3/2 becomes 



5^3/4 



«3/2 < 1-9 X 10-^ a, 



Ms 



IQio GeV 



-3/8 



V8 /Y- 



IQio GeV J \ 10-16 



BBN\ 1/2 



3/2 



(50) 



and the one on the branching ratio 
B3/2 < 9.8 X 10-12 ^- 



M5 



-3/4 



Mx 



1010 GeV 



,/4 ^ YjJBN - 
10-16 



(51) 



, 1010 GeV 

As in SC [cf. Eqs. ([13|) and (jUj) ]. a very low value for Q!3/2 (or equivalently for S3/2) is also required 
in the high-energy regime. This implies that the decay of the heavy scalar into gravitinos is very 
constrained. We will comment more on this in the conclusion. 

In the following we will discuss in more detail the additional cosmological constraints on the 
gravitino amount produced at the reheating by the X decay. Namely, we will use BBN constraints 
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on the amount of unstable gravitinos to impose bounds on the relevant quantities. Afterwards, we 
will discuss the case when the gravitino is stable and use the WMAP limits on the amount of cold 
dark matter to derive constraints. 



B. Unstable gravitino and BBN constraints 

Let us start with the unstable gravitino case. In conventional scenarios, the gravitino mass is 
expected to be comparable to the masses of the supersymmetric partners of the standard model 
particles and, therefore, m3/2 ^ a few TeV in order to solve the gauge hierarchy problem. Since 
the gravitino coupling to matter is suppressed by the Planck mass Mp , its lifetime is ~ 
Mp/my2 ~ 10^(100 GeV /m^^2)^ s. For a gravitino mass > 100 GeV, it decays during or soon 
after the BBN epoch. Therefore, if gravitinos are overabundant, the high energetic photons and 



hadrons emitted in their decays may destroy t. 
BBN. However, it has been recently pointed out 



le light elements and hence spoil the success of 



331 ] that the hadronic decay gives a more stringent 



constraint on the abundance of gravitinos, because mesons and nucleons produced in the decay 



and subsequent hadronizatioii 
that come from BBN are 



arocesses significantly affect BBN. For rn-3/2 = 1 TeV, the bounds 



10| 



4 X 10-1^ for Bh = 1 (Case I) , 

5 X 10-^'^ for Bh = 10-3 



Y,T<{ ^ T^^'^ (52) 



where Bh denotes the hadronic branching ratio of the gravitino. While, for = 10 TeV one 



has 



lOfi 



2 X 10-^^ for Bh = l , 
2 X IQ-^^ for Bh = 10-3 (Case II) . 



The region ^ 10^ TeV is disfavored since large gravitinos masses do not solve the gauge 

hierarchy problem, as mentioned before. The BBN bound disappears for m^/2 ~ 10^ ~ 10^ TeV 
and, for this mass range, the constraint comes from the abundance of the relic LSP produced by 
the gravitino decay ^J^. 

In Fig. [3] we present the contours of constant 1^3/2 = ^3/^^ ^ov several values of M5. We 
present the contours for 171-^/2 — 1 TeV and Bh = 1 (case I), for which ^^^^ ^ 4 x 10-^^, and 
m3/2 = 10 TeV and Bh = 10"^ (case II), for which Igyf^ < 2 x lO^^^^ the right panels we 
show two sets of curves, one for at = 1 and the other for at = W^^. Once more we have fixed 
rrig = 1 TeV in the computation of the thermal contribution to the gravitino abundance. Some of 
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a^^^ (GeV) (GeV) 



FIG. 3: (Color online) Contours of constant Fg/s = Yf^^^ for Afg = 10^ 10^ 10^°, lO^i, lO^^, IQi^ GeV 
(from lighter to darker green). In the top panels we have Y^^^ = 4 x 10^^^ (^3/2 = 1 TeV and Bh = 1), 
and in the bottom ones Y^^^ — 2 x 10^^^ ("t-3/2 = 10 TeV and Bh = 10^^). In the right panels we show 
two sets of curves for which we fixed at = I (full lines) and at = 10^^ (dashed lines). For 7713/2 — 1 TeV, 
the A/5 ~ 10^^ GeV curves already reproduce the results from SC, while for 7713/2 = 10 TeV we obtain the 
SC results only for M5 - 10^^ GeV. 

the curves are broken due to the upper limits on Trh and Mx [cf. Eqs. (f2T|) and (f22|) . respectively] 
which come from 7723/2 > H. 

One should notice that for 7713/2 = 1 TeV, the curves corresponding to M5 > 10^^ GeV already 
reproduce the results from SC, while for 772,3/2 = 10 TeV we obtain the SC results only for M5 > 
10^^ GeV. The value of M5 below which the BC effects affect the production of gravitinos can be 
easily obtained from Eq. ()47p . One has 



M5 < 3.3 X 10^" GeV , for 7773/2 = 1 TeV , 



(54) 
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and 



Ms < 9.7 X 10^^ GeV , for mg/a = 10 TeV . (55) 

For 5D Planck masses above these values, the upper bounds on Tj-^, Mx and 03/2 derived for 
SC are valid. For the reheating temperature one has [cf. Fig. [3] (left) and Eq. (j41|) ] 

' 4.7 X 10^ GeV (I) , 
2.1xl09GeV (II), 
while for the heavy scalar mass one finds [cf. Fig. [3] (right) and Eq. (j42p ] 



Trh< 



(56) 



Mx< 



4.1 X lO^Oa^^/^GeV (I), 

1.2 X IQi^a^^/^ GeV (II) , 



(57) 



and for 0^/2 [cf. Fig. [3] (right) and Eq. ([l3] 



"3/2 ^ 



1.3 X lO-^a^^ 



Mx 



-1/4 



9.4 X 10"* a 



4 ^1/2 



IQio GeV 

Mx 



(I) , 
(11) 



(58) 



IQio GeV^ 

The lower bound on the reheating temperature presented in Eq. (j49p reads in this case 



3.4 X 10^ GeV ( (I) , 



Trh > <! 



105 GeV J 



6.8 GeV 



03/2 MxV 



(59) 



,10^; 

For lower values of the 5D Planck mass, the disappearance of the upper bounds on Trh and 
Mx coming from BBN limits on the amount of gravitinos will allow the Universe to reheat at a 
larger temperature and, at the same time, the scalar to have a larger mass and a larger decay rate 
into gravitinos than the ones allowed in SC. The upper bound in AIx is simply the one given in 
Eq. (I22p . while the reheating temperature, in addition to T^h > Tt, should satisfy 



2/3 



6.1 X 10^ GeV f^^) 

03/2 Mx\ 



3.5 X 10^ GeV 



lo^Gevy Vio^°Gevy < 



105 GeV 



4.2 X 10^ GeV 



Ms \ 
1010 GeV J 
Afs A 



3/4 



3/4 



(I), 



(60) 

obtained from Eqs. (08]) and (fT9]) . respectively. In this regime, the upper bound on the coupling 
constant will grow with Mx, in contrast with the SC scenario. From Eq. (jSOp we obtain 



"3/2 ^ 



5 ^3/4 



1.2 X 10-5 a 



^.4 X 10-4 ^ 



M 



X 



1010 GeV 

Mx 
1010 GeV 



1/8 



Afs 



1010 GeV 

Ms 
1010 GeV 



-3/8 



-3/8 



(I) , 

(II) . 



(61) 
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Before concluding this section, it is worth commenting on the possible additional constraints 
coming from the LSP abundance, which is stable if R-parity is conserved. In this case one must 
impose that the abundance of LSP does not overdose the universe, as it contributes to the dark 
matter content. The LSP can be produced both from gravitino and heavy scalar decays. The 
LSP abundance from gravitino decays does not give additional bounds provided that 777-3/2 = 
100 GeV — 10 TeV. The constraints coming from the production of the LSP from heavy scalar 
decays will be model dependent. In Ref. the bounds arising from considering the neutral wino W 
as the LSP, which maximizes the annihilation cross section and hence leads to the most conservative 
bounds, were derived in SC. In order to avoid the overclosure of the universe, a lower bound in Mx 
is found, namely Mx ^ 2 x 10^ GeV for Af^ = 100 GeV, which becomes more stringent for larger 
wino masses. Since the brane effects are only present for heavy scalar masses Mx ^ C(^ M5 [cf. 
Eq. (f22l) ]. one should expect that for sufficiently high values of M5 these bounds will not change. 



C. Stable gravitino and dark matter 

In the case of a stable gravitino, when the gravitino is itself the LSP with exact R-parity, one 
should check whether the contribution of gravitinos to the energy density of the universe does not 
exceed the observed matter density limit. From the gravitino abundance, we can estimate their 
contribution to the closure density 

^3/2h^ = ^3/2 Y3/2 So h^Pc^ ■ (62) 

Here pc = 3iJ^M|,/87r = 8.07 x IQ-^'^h'^ GeV^ is the critical density and sq = 2.22 x lO^^s GeV^ 
is the present entropy density. The total fraction of gravitinos contributing to dark matter, ^^3/2, 
is the sum of the two contributions 

^^3/2 = ^3/2 + ^3/2 ' (63) 



the first coming from the abundance of gravitinos due to the decay of the heavy scalar, Eq. ()32p . 
and the second from thermal production, Eq. (p6|) . To avoid overclosure one must require 



^3/2h^ < ^DMh? , (64) 
where we will use the WMAP bound on the matter density of the universe jsj], 

= 0.143 . (65) 
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If ^^3/2 = ^DM, the gravitinos constitute (all) the dark matter of the universe. 

Let us now consider the two contributions to the gravitino dark matter. From the decay of the 
heavy scalar one has 



n^/2h^ = 2.8 X 10^^ 



= 4.1 X 10"' 
and from the thermal production 




Trh \ 

106 GeV J 



nfi^h^ = 2.8 X 10"^ 



m 



3/2 



10-1 GeV 



In the SC case one can write 



''3/2"- 



2.3 X 10"^ at 



3/2 



m 



3/2 



^ifl^h^ = 1.6 X 10-2 



10- V V 10-1 GeV 

A{Trh; 171^/2, nig 



10 



-14 



GeV 



' vth - 
-^3/2 

10-9 



Mx 
lOio GeV 

"^3/2 



1/2 



10-1 GeV 



Mx \ 
1010 GeV J 



3/2 



while in the BC regime one has 



-2^1/2 



'3/2' 



2.8 X 10"^ a 



4.4 X 10^2 



B 



3/2 



m3/2 



10-4 y V 10-1 GeV 
A{Tt; 1713/2, mg 



Mx 



1010 GeV 



-1/4 



3/4 



10 



-14 



GeV" 



lOio GeV J 

^3/2 \ f Ms y/' 

10-1 GeV/ VlO^^GeVy 



(66) 

(67) 

(68) 
(69) 

(70) 
(71) 



Notice that, in both BC and SC, il^g^^ proportional to 777,3/2, while 0.^2^^"^ is inversely pro- 

This means that for large gravitino masses the 



~ m 



3/2- 



portional to 777.3/2, since ^(T^; 7773/2, 777^) 
X decay contribution will dominate, while for small values of 7773/2, the main contribution comes 
from thermal production. On the other hand, the dependence of 0^2^^ Mx is different in BC 
from that of SC. Once more, this implies that in the brane high-energy regime, there is no upper 
limit on the heavy scalar mass coming from the gravitino dark matter fraction. 

One should notice that the momentum of the gravitinos produced by X decay can be much 
larger than its mass since P3/2 — Mx/2. Therefore, they can behave as warm dark matter and 
influence the structure formation: with a non negligible velocity, and before the time of matter- 
radiation equality, they can freely stream out of over dense regions and into underdense regions, 
and this can erase the small structures that are observed today. This leads to an upper bound on 
the present dispersion velocity of the gravitino v'^^^ 



35l . l36l . l37l | . The present velocity of dispersion 



of the gravitinos produced by the X decay is estimated to be 



35| 



7,0 
^3/2 



1 Mx To ( g,{T^ 



2 7773/2 Trh \g*{Trh 



1/3 



(72) 
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where Tq = 2.35 x 10 GeV is the present photon temperature and 5*(To) = 43/11. Bounds on 



this velocity are obtained by means of the power spectrum inferred from 



wnan-a together with 



39|. From Ref. 



• 3 



cosmic microwave background radiation and galaxy clustering constraints 
it was found [9] that v'^^^ < 4 x 10"^. 

From the limit on the present velocity of dispersion one can estimate an approximate upper 
bound on the branching ratio of the decay channel into gravitinos i?3/2- Considering the contribu- 
tion from the heavy scalar decay to the amount of gravitinos we can write 

nf/. = 5%2^»0P.-', (73) 



from which 



-3/2 = 1-26 X 10-^^ . (74) 



Assuming n^jr^h? = Q.^j^.^h? , one gets the approximate bound 



i33/2 < 4.5 X 10-^ (75) 



Notice that this bound is independent of the background cosmology. One can also compute the 
velocity of dispersion for gravitinos produced thermally, however one finds that for 711,3/2 ^ 100 keV 
they have a neglig we will consider masses well above this 

bound. The above constraint in -B3/2 was determined assuming that all dark matter is constituted 
by gravitinos, ^^^9 = ^dm- However, if gravitinos are only a fraction of dark matter, then the 
bound becomes weaker, and can even be evaded if 

i7f/2 <0.12J^BM, (76) 



as derived in Ref. 



In Fig. m we present the contours of constant 1^3/2/1^ = ^DAih'^, in the Mx — plane, with 
^DAih^ given by the WMAP bound in Eq. (j65p . The different (green) full lines correspond to 
various values of -63/21 where we fixed at = 1 and rrig = 1 TeV. The allowed range of masses 
lies to the left and/or below the contour lines. Moreover the gravitino dark matter is only viable 
to the right/above the line which corresponds to the upper bound on the branching ratio shown 
in Eq. ([75]) . -B3/2 = 4.5 x 10~^. The yellow shaded area is excluded because of the condition 
H{Trh) < 1^3/2 [cf- Eq. (f22l) ]. The horizontal line represents the value of Mx = (2/a^) M5, for 
which the transition from the high to the low energy regimes occurs. The contours are shown for 
three different values of the 5D Planck mass, namely, M5 = 10^, 10^^, 10^'^ GeV. 
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lO"" 10^ 10° 10^ 10"' lO"'^ 10° 10^ 10 ' 10^ 10° lo'' 

m,,, (GeV) m (GeV) m (GeV) 



FIG. 4: (Color online) Contour of constant ^13/2^1^ = ^DAih'^ = 0.143 for B3/2 = 10"^, 4.5 x 
10"'*, 10"**, 10"^ (darker to lighter, or left to right, green full lines) in the Mx — m-3/2 plane. The dif- 
ferent panels correspond to various values of M^. We have chosen rrig = 1 TeV and at — 1- Also shown 
is the contour for 0.^/2^^ = Q-^'^.VLjjMh^ (blue dashed line), when i?3/2 = 10"^. The yellow shaded area is 
excluded because of the condition H{Trh) < [cf. Eq. (|22p ]. The horizontal line represents the value of 

Mx for which the transition from the high to the low energy regime occurs. 

In SC there is an upper limit on the mass of the heavy scalar Mx coming from the bounds on 
the gravitino abundance. The maximum value of Mx is obtained when both contributions to the 
gravitino production, from the heavy scalar decay and thermal scatterings, are equal; hence, at 
that point, one can write that the total abundance is $^3/2 = '^\J^^/2 ^3/2 ^ ^dm and obtain 

Mx < 7 X 10^ GeV a~l^ . (77) 

We have taken into consideration that, for the gravitino mass in the range 10~^ — 10^ GeV , 
A{Trh'-,'ni^l2-,inng) ~ IG^^'^m^J^ is a good approximation. As said before, in the high energy regime 
there is no such a bound, and the allowed values for Mx are only constrained by the condition 
H{Tj.h) < fn^i2 [cf. Eq. ([22]) ]. We see that, for low values of M5 and a fixed it is possible to 

have larger values of both Mx and m3/2 and yet satisfy the bounds from dark matter. For example, 
for M5 ~ 10^ GeV and sufficiently low branching ratio, it is possible to have Mx ~ lO^'^ GeV while 
having m3/2 ~ 100 GeV. 

If one now takes T^h as the free parameter one can derive the contours of constant $73/2^1^ = 
^D&ih'^ in the Tj-h — ?^3/2) shown in Fig. [5] for different values of M5. Here the different (green) full 
lines correspond to various values of 03/2^3/2) and the gravitino dark matter is only viable to the 
right/above the line which corresponds to the upper bound on the branching ratio -B3/2 ^ 4.5 x 10~^ 
[since 03/2 ^ C(l)]- Once again one can distinguish the two behaviors for SC and high-energy BC 
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10"' 10" 
m,„ (GeV) 



FIG. 5: (Color online) Same as in Fig. [H but in the Trh ~ '7^3/2 plane and for 013/253/2 = 10 4.5 
10""*, 10~^, 10"'' (darker to lighter, or left to right, green full hues). 



regimes. The contributions to the gravitino dark matter read as 



n^,,h' = 1.9 X 10-2 ( 



2/3 



"^3/2 



10-1 GeV 



1/3 



nf^^h'^ = 2.8 X 10"^ 



-^{Tt; 1713/2, m-g] 
10-14 GeV-i 

in SC hmit. In the BC regime we obtain 



m 



■3/2 



10-1 GeV 



106 GeV^ 

Trh \ 

106 GeV J 



^3/2^^ 



3.5 X 10" 



«3/2 -63/2^ 



nfj^h^ = 5.5 X 10"^ 

As expected, the bound 



10-6 
^(Tj; 7713/2, mg 



10 



-14 



GeV- 



"^3/2 \ / Trh \ 

10-1 GeV/ V 106 GeV y 
"^3/2 \ / Tt 
10-1 GeV/ V 106 GeV 



-1/3 



1010 GeV 



(78) 
(79) 

(80) 
(81) 



Trh < 1.8 X 10^ GeV 



1 /2 

"3/2-^3/2"^ 



10 



-5 



(82) 



derived for the SC case in order to avoid the overclosure of the universe by dark matter disappears 
when we consider sufficiently low 5D Planck scales M5. 

If -B3/2 > 4.5 X IQ-^, then the gravitinos can only contribute to a fraction of the total amount of 
dark matter. In Figs. [4] and [5] the blue-dashed line represents the contour for 1^3/2/1^ = 0.12ri£)A//i^, 
the maximum contribution of gravitinos (warm) dark matter that does not spoil the structure 



formation, for B 



3/2 



10" 



We can see that in this case the bound becomes more stringent. 



Finally, in Fig.[6]we show the contours of constant 0,^/2^^'^ = ^Dhih'^i in the T^/j — a3/2Mx plane, 
for two different gravitino masses: 7713/2 = 1, 0.1 GeV (full and dotted bold lines, respectively) 



and M5 = 10 , 10 , 10 GeV. Also shown are the contours of constant -B3/2 (dash-dotted lines) 
for a fixed 0:3/2 = 1. It is clearly seen from the figure that, in the BC regime and for sufficiently 
low values of M5, higher reheating temperatures are allowed for large Mx- 
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"3/2 '^''x (°''^) °-3l2 l^"^' "3/2 '^''x (°''^' 

FIG. 6: (Color online) Contours of constant Vlj,/2h?' ~ flnAih^ ~ 0.143 for m^^2 ~ 1 GeV (full bold line) and 
"13/2 — 0-1 GeV (dotted bold line) in the Tr/i — a3/2^^x plane. The dash-dotted lines represent the contours 
of constant -B3/2 = 1, 10^^, 10^^", 10^^^, 10^^" (from below to top), and the full line is the contour for 
^3/2 — 4.5 X 10^'*. In the computation of the branching ratio we fixed 0:3/2 — 1- The different panels 
correspond to various values for M5. 

IV. DISCUSSION AND CONCLUSION 

We have studied the gravitino production in the braneworld cosmology context. Our framework 
is the RSII construction in which the Friedmann expansion law is modified with a quadratic term in 
the energy density. As an approximation, we considered the gravitino as a field localized on the UV 
brane and used 4D supersymmetry and the usual Boltzmann equation to describe its production. 

In Figs. [7] and [8] we show the points, in the T^h — Mx, T^-h ~ -^5 and at — 03/2 planes, that 
satisfy the constraints discussed in detail in the other sections. We have allowed a variation of the 
couplings and masses in the ranges 03/2 < 1, 10 ^ < a* < 1, 10^ GeV < Mx < 10^^ GeV 

and 10^ GeV < M5 < 10^*^ GeV. For the stable gravitino case one has 10""^ GeV < 7713/2 < 10^ GeV. 
The light and dark green points correspond to the SC and BC high energy regimes, respectively. 

For the unstable gravitino case, one can derive the range of M5 for which the brane correction 
to the standard expansion rate is relevant. The upper bound for this range is dependent of the 
gravitino mass and we have obtained M5 < 3.3 x 10^'^ GeV and M5 < 9.7 x 10^^ GeV, for 
m3/2 = 1, 10 TeV, respectively. 

In the BC regime the upper bounds on the mass of the heavy scalar and the reheating tem- 
perature, coming from the gravitino abundance, disappear. Thus, with the inclusion of the brane 
correction, the parameter space with the correct gravitino abundance is enlarged. Both, larger 
reheating temperatures and larger masses for the heavy scalar are allowed. This is more evident 
for a smaller gravitino mass and a larger hadronic branching ratio. Although Trh and Mx are not 
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FIG. 7: (Color online) Parameter space allowed by the gravitino abundance BBN bounds, in the T^h — Mx, 
Trh — and at — 03/2 planes. In the top panels we have Y^^!^^ = 4 x 10^^^ ('713/2 — 1 TeV and Bh = 1), 
and in the bottom one Y^^^^ = 2 x 10"" (7713/2 = 10 TeV and Bu = 10"^). The light and dark green 
points correspond to the SC and BC high energy regimes, respectively. 



constrained by the gravitino abundance limits, they are bounded by the fact that H > 7713/2 for 
the gravitino decay to be possible. 

On the other hand, in BC, like in SC, one also needs a small coupling 03/2 and/or a small 
branching ratio -63/2- However, for a fixed at, it is possible to have a slightly larger 03/2 in BC. 
The required lower value for 03/2 does not allow to solve the moduli problem by simply increasing 
the moduli mass, since, for the range of gravitino masses considered, Q3/2 = 0{1) is excluded. It 
also poses severe constraints to inflationary models. In order to achieve a small a3/2 one needs a 
VEV of X (the inflaton in this case), much smaller than Mp; on the other hand, one can decrease 
B-^I2 if one increases at, by demanding interactions of X with the other MSSM particles with 
strength much larger than Mp^ . A lower value for a3/2 can be achieved in some new and hybrid 
inflationary models, and in the simplest chaotic inflationary models, for which Gx vanishes in the 
vacuum (see Refs. [a, 0, 1^) one obtains Q3/2 = and then no gravitinos are produced in the decay 
of X. 

In the stable gravitino case, similar bounds can be derived. Notice however that the overclosure 
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FIG. 8: (Color online) Parameter space allowed by the dark matter bounds, $13/2^^ < ^DAih^ = 0.143, 
and ^3/2 ^ 4.5 x 10^* , in the Trh ~ Mx, Trh — M5 and at — a3/2planes. The gravitino mass runs in the 
interval 10^^ — 10'^ GeV. The light and dark green points correspond to the SC and BC high energy regimes, 
respectively. 

constraint alone allows to have -B3/2 ~ 10~^. Therefore, it would be possible to solve the moduli 
problem. But, if one takes into account the warm dark matter constraint, then one needs -B3/2 ^ 
10~^, which means that Q3/2 ^ 0.1 at (cf. Fig. [8|). In Ref. its was suggested that such a 
suppression in the coupling might be possible in some moduli stabilization mechanism. In that 
case it would be possible to have gravitino warm dark matter from the decay of a moduli field with 
mass as large as 0(10^^) - 0(10^^) Qgy, for a gravitino with mg/g < 3 GeV and M5 ~ 0(10^) (cf. 
Fig. S]). In SC, the upper bound for the mass of the moduli would be one order of magnitude bellow 
the BC one. Moreover, we should point out that the warm dark matter constraint disappears if 
the gravitino contributes with less than 12% to the total dark matter density. 

Finally, it would be interesting to further investigate the gravitino production in specific brane 
inflationary models. The gravitino production analysis presented here could also be extended to 
include the KK gravitino production process. 
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